Transient effects and reconstruction of the energy spectra in 
the time evolution of transmitted Gaussian wave packets 
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We derive an exact analytical solution to the time-dependent Schrodinger equation for trans- 
mission of a Gaussian wave packet through an arbitrary potential of finite range. We consider 
the situation where the initial Gaussian wave packet is sufficiently broad in momentum space to 
guarantee that the resonance structure of the system is included in the dynamical description. We 
demonstrate that the transmitted wave packet exhibits a transient behavior which at very large dis- 
tances and long times may be written as the free evolving Gaussian wave packet solution times the 
transmission amplitude of the system and hence it reproduces the resonance spectra of the system. 
This is a novel result that predicts the ultimate fate of the transmitted Gaussian wave packet. We 
also prove that at a fixed distance and very long times the solution goes as t~^^'^ which extends 
to arbitrary finite range potentials previous analysis on this issue. Our results are exemplified for 
single and multibarrier systems. 
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I. INTRODUCTION 

The transmission of wave packets through one- 
dimensional potentials is a model that has been of great 
relevance both from a pedagogical point of view, as dis- 
cussed in many quantum mechanics textbooks, and in 
research, particularly since the advent of artificial semi- 
conductor quantum structures [uja. There are studies 
on the dynamics of tunneling ^UHj H 0, B or work 
on controversial issues, as the tunneling time problem 
[Tol . [ill. IT^ ll^ . and on related topics as the Hartman 
effectjlllljor the delay time [11 [13, El- Most time- 



dependent numerical studies consider Gaussian wave 
packets as initial states [1, [1, [l^ [l^l , though in some re- 
cent work, the formation of a quasistationary state in the 
scattering of wave packets on finite one-dimensional pe- 
riodic structures involves also some analytical considera- 
tions [2]| . Analytical approaches have been mainly con- 
cerned with cutoff quasi-monochromatic initial states in a 
quantum shutter setup [1, [l^, [13, [H, [l^l . In recent work, 
however, analytical solutions to the time-dependent wave 
function have been discussed using initial Gaussian wave 
packets for square barriers [^[Sll, delta potentials [6l[25| 
and resonant tunneling systems near a single resonance 
i- 

We obtain an exact analytical solution to the time- 
dependent Schrodinger equation for transmission of an 
initial Gaussian wave packet through an arbitrary po- 
tential of finite range. We refer to the physically relevant 
case where the initial Gaussian wave packet is sufficiently 
far from the interaction region so that the correspond- 
ing tail near that region is very small and hence may 
be neglected. Since the infinite limit of a very broad 



cutoff Gaussian wave packet in configuration space, i.e, 
that leading to a cutoffplane wave, has been discussed 
analytically elsewhere 3], we focus the discussion here 
to cases where the initial cutoff Gaussian wave packet 
is sufficiently broad in momentum space so that all the 
resonances of the quantum system are included in the dy- 
namical description. We demonstrate that the profile of 
the transmitted wave packet exhibits a transient behav- 
ior which at very large distances and long times may be 
expressed as the free evolving wave packet modulated by 
the transmission amplitude of the system. To the best of 
our knowledge this is a novel and interesting result. We 
also analyze the transmitted solution at a fixed distance 
away from the potential at very long times, and find that 
it behaves as i^"^/^. Our result generalizes to arbitrary 
potentials of finite range previous analysis involvin g sp e- 
cific potentials models and numerical calculations [26|. 
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The paper is organized as follows. In section II, some 
relevant aspects of the formalism of resonant states are 
reviewed and a formal solution for the transmitted pulse 
is given as an expansion in terms of these states and the 
corresponding complex resonance poles. In Section III 
the analytical expression for the transmitted Gaussian 
wave packet is derived and some limits are discussed. 
Section IV refers to some examples, specifically single 
barrier and double and quadruple barrier resonant tun- 
neling systems are considered and a subsection provides 
some remarks concerning the tunneling time problem. 
Section V gives the concluding remarks and, finally, the 
Appendixes discuss, respectively, the analysis of the ef- 
fect of the cutoff in the solution, and a general method 
to calculate the complex poles of the transmission ampli- 
tude. 
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II. RESONANCE EXPANSION OF THE 
TIME-DEPENDENT SOLUTION 

Let us consider the time evolution of an initial state 
ii){x, 0) of a particle of energy Eq = h^kl/2m, approach- 
ing from X < a. potential V{x) that extends along the 
interval < x < L. The time-dependent solution along 
the transmitted region x> L reads |22, i27i | 



V'(x,t) 



1 



/2tt J- 



dk (/)o(fc)t(fc)e 



ikx—ihk^t/2m 



(1) 



where t(fc) is the transmission amplitude of the problem 
and (/>o(fc) is the Fourier transform of the initial function 
i:{x,0). 

One may write the transmission amplitude in terms of 
the outgoing Green function G^{x, x'; k) of the problem 
as 



t(fc) = 2ikG+{0,L;ky 



-ikL 



(2) 



It is well known, that the function G^{x,x';k), and 
hence the transmission amplitude t(fc), possesses an infi- 
nite number of complex poles k„, in general simple, dis- 
tributed on the complex k plane in a well known manner 
[1^. Purely positive and negative imaginary poles k„ = 
i"/n correspond, respectively, to bound and antibound 
(virtual) states, whereas complex poles are distributed 
along the lower half of the k plane. We denote the com- 
plex poles on the fourth quadrant by k„ = a„ — i/3„. It 
follows from time reversal considerations [2^ that those 
on the third quadrant, fulfill K_n = — The 



complex poles may be calculated by using iterative tech- 
niques as the Newton- Raphson method [33| , as discussed 
in the Appendix. Usually one may obtain a resonance 
expansion for t(fc) by expanding G^{0, L; k) in terms of 
its complex poles and residues [3Q| . Here we find more 
convenient to expand instead 0^(0, L; fc) exp(— i/cL) to 
obtain, 



t(fc) = ^fc y 



(3) 



where the residues r„ are given by 

_ ■u„(0)u„(i) 



(4) 



The functions Un(x) appearing in Eq. dH) satisfy the 
Schrodinger equation to the problem with complex eigen- 
values En — ?i^k^/2to —£n — and obey the purely 
outgoing boundary conditions 



dx 



Unix) 



x=0 



dx 



Unix) 



iK„UniO) 

normalized according to the condition 

uliO) + uliL) 



u^ix)dx + i- 



2hir, 



iK„u„iL) 
(5) 

(6) 



Substitution of Eq. ^ into Eq. H]) yields, 



7 °° r°° k 
iPix,t)^^ y r„e-'«"^ / dk- 0o(fc)e*'="-''^'= 'Z^"* (7) 

V 27r J-oo ~ ^n 



I 

It is convenient to make use of the identity /c/(fc— k„) = and ipnix, t) is given by 
1 + f^n/ik — Kn) to rewrite the solution given by Eq. ([7]) 

as '^nix^t) — ZT^k^G ^ X 

i;ix,t)=C^P^ix,t)+ y i^nix,t), (8) r J^M^^^kx-^hkh/2m^ 

nt^oo J-^V2^k- Kn 

where C is a constant that depends only on the potential 

through the values of the r„'s and k„'s. III. ANALYTICAL SOLUTION FOR A CUTOFF 

GAUSSIAN PULSE 

Consider now a particle described initially by a cutoff 
Gaussian wave packet 

( Aoe-(^-"^<^)'/'*'"'e*'=o^, a; < 
^fix,t) = -1= jykMk)e"'''/'-, (10) ^ I ^ ^ ^ ' 



C = z y r„e-''^"^, (9) 

n— — oo 

tp^{x^t) stands for the free wave packet solution 
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where Aq is the normaHzation constant, Xc, cr and ko 
are, respectively, the center, the effective width and the 
wavenumber corresponding to the incident energy Eq of 
the wave packet. 

In order to calculate the free evolving wave packet 
given by Eq. ()10|) and the integral term on the right- 
hand side of Eq. one needs to know the Fourier 
transform of the Gaussian cutoff wave packet. This is 
given by [2^ . 



4>o{k) = Aouj{iz), 



where 



An- 



/27r 



where z and zq are given by 



and 



Xq 

2^ 



zo 



i{k — fco)cr. 



V2a 



(13) 



(14) 



(15) 



(16) 



and U!(iz) is the Faddeyeva function [3lLl32l|. 

Let us place the initial wave packet along the region 
Xc < 0. As pointed out above, here we shall be concerned 
with the physically relevant situation where the tail of 
the initial Gaussian wave packet is very small near the 
interaction region. It is then convenient to consider the 
symmetry relationship of the Faddeyeva function [3ll.[3^. 



Lo(iz) = 2e^ — uj{—iz) 



(17) 



and follow an argument given be Villavicencio et.al. for 
the free and 5 potential cases [1^. These authors obtain 
that provided 



2^ 



> 1; Xc< 0, 
one may approximate Lj{iz) as 

u){iz) ~ 2e^ . 



(18) 



(19) 



In Appendix A we show that the above approximation 
holds also for the general case of finite range potentials. 
Using Eqs. and HH) into Eq. HU]) leads to an an- 

alytical expression for the free evolving cutoff Gaussian 
wave packet [1^, that we denote by ipj{x, t), that is iden- 
tical to the exact analytical expression for an extended 
initial gaussian wave packet [25| , 



1 



-j^ ^{{kox — hk^t /2m) 



(2^)1/4 ^1/2 ^1 + it/r 

{hko/m)tY 



exp 



X Xr 



40-2 [I + it/r] 



(20) 



where 



2ma 



(21) 



Let us now substitute Eq. ([T9)l into the integral term 
in Eq. (fTT|) to obtain 



. dk 
2Ao 



^ikx — ihk^t/2m ^22) 



, V 27r k K, 

Feeding the expression for z appearing in Eq. (jl5p into 
Eq. dm allows to write Tp!;^{x,t) as 

ij^{x,t) = irnKne-"^-^ X 

De^k„^-^hklt/2mj^j^^y,^^^ (23) 



where 



D - -2i(27r)i/'*V^> 



(24) 



that follows using ^Jijj{izQ) = exp(zQ/2) erfc(zo) with 

v/erfc(zo) ~ 2 for zq < -1 [sil, [13], and M{y'^) stands 
for the Moshinsky function, defined as [1, [2^ , 



. oo ikx' —ihk^t' /2m 

M[y'^) ^ ^ I dk 



= ie™-"/2^*'c.(zy;). 



(25) 



with 



hko 



x' — X — Xc 1, t' — t ~ ir, k' — Kn — ko, (26) 

TO 



and the argument y'„ of the Faddeyeva function uj{iy'„) 
reads, 



t/4 



2nt' 



TO 



(27) 



The Moshinsky function is usually calculated via the 
Faddeyeva functions for which well developed computa- 
tional routines are available [l^j , Substitution of Eq. (1^5)) 
into Eq. ([5]), allows to write the time-dependent trans- 
mitted solution as, 
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OO 

ij''{x,t) =Cij'}(x,t)+iDe'^°''^'^''o'/^"' rnKne-"'-^M{y'J. (28) 

n— — OO 



Notice, in view of the definitions for x' and t' in Eq. 
(PS)) and of T given by Eq. (|2ip . that the argument of 
the second exponential term on the right-hand side of 
Eq. may be written as —imx'^/2'ht'. This allows to 
write Eq. ^ for ipl{x,t) as 

1 1 ^i{kQX — hkQt/2m)^ — imx''^/2ht' 

(29) 



and Eq. (pS)) . using Eq. (pS)) . alternatively, in the more 
convenient form as 



i;^{x,t)=i,'}ix,t) 



E 



(30) 



One should emphasize that Eqs. ([29]) and (|30p hold pro- 
vided the condition given by Eq. (|18p is satisfied. 

Clearly, as a consequence of the approximation given 
by Eq. (fTO|) . the solution ip'^{x, t) does not vanish exactly 
as t — > 0. There remains a small value proportional to 
the tail of the free solution. 



A. Long-time behavior of ip'^{x,t) 

Let us now analyze Eq. (ISOp at asymptotically long 
times, i.e. much larger than lifetime r of the system, 
for a fixed value of the distance x — x^. In such a case, 
one sees from Eq. (|27p . that the argument of the 
Faddeyeva functions behaves as 

« _e— /4 ^^^^^ (31) 
V 2m 

and hence becomes very large as time increases. For 
proper poles k„ = a„ — i.e., q;„ > /3„, the Faddeyeva 



function behaves as [3l|, [13] , 

«;(zy;)«e^" +4? (32) 

Using Eq. ([5T|) it follows that the term exp(j/^^) vanishes 
exponentially with time. On the other hand, for poles 
K-n = —K^, seated on the third quadrant of the k plane, 
the Faddeyeva function behaves in a purely nonexponen- 
tial fashion as on the right hand-side of Eq. ([5^ [sil . [s^ , 
namely 

«;(z2/U)«4if^-4l^V-"' (23^ 

T^^'' \y-n y'tnj 

One sees therefore, that for sufficiently long times the 
full set of resonance poles behaves nonexponentially. Us- 
ing Eqs. dSl), dSll) and §3^, one may write Eq. dSO]) at 
asymptotically long times as, 



E 



E 



-iKjiL 



3/r 
y'l. 



(34) 



Substitution of Eq. (|3T|) into Eq. (p4)) . one sees that the first one and hence one obtains that ■ip'^{x,t) behaves as, 
second term on the right-hand side cancels exactly the 

ip^ix^t) -jj^; X ^ Xd, t'y^T. (35) 
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It follows from the above expression that the correspond- 
ing probability density goes as This long-time 
behavior of the probability density as an inverse cubic 
power of time has also been obtained with other poten- 
tial models and initial states, including numerical calcu- 
lations of Gaussian wave packets colliding with square 
barriers [1^ . As pointed out in Ref . [1^ the above long- 
time behavior for the probability density is consistent 
with the definition of the dwell time as a physical mean- 
ingful quantity. 



B. Asymptotic behavior of ^l}°'{x,t)/xp'f{x,t) 

There is another asymptotic limit involving the trans- 
mitted wave packet solution given by Eq. pO]) . This 
refers to the hmit of tp"-{x,t)/tlj'j;{x,t) as a; ^ oo and 
t — > oo. Previous analysis regarding the time evolution of 
forerunners involving cutoff initial plane waves show that 
at very large distances and long times, x/f — > (hk/m) 
[36| . This suggest a similar behavior for the transmit- 
ted Gaussian pulse. Hence as x and t attain very large 
values, one may write the argument y'„ of the Faddeyeva 
function, given by Eq. ([27]) as 

y'n^e''^^\h^[k-Hn]t'^^; X={hk/m)t, t^OO, 

V zm 

(36) 

where the relationships given by Eq. ()26p have been 
used. It follows then, using the leading 1/y' terms 
in Eqs. (132]) and ((33)) . that at very long times the 
term ^/l + it/T Lo{iy'^) appearing in Eq. ([50)1 tends to 
i/[{TT^/'^a){k — K„)]. As a consequence, one may rewrite 
Eq. (|30)) , for very large values of x and t as 

^!(^=t(fc); x^{hk/m% t»T (37) 

iPj[X, t) 

where Eqs. ([9]) and ([3]) have been used. Equation 
P7p provides an analytical demonstration that at very 
large distances and times, ^'^ {x , t) / tpj {x , t) reproduces 
the transmission amplitude of the system, and hence 
\ip°-{x,t)/->p'j{x,t)\'^ = T{E) vs E, the corresponding 
transmission energy spectra of the system. 



IV. EXAMPLES AND DISCUSSION 

In order to exemplify our findings, we consider three 
tunneling systems involving typicalparameters of semi- 
conductor AlxGai^xAs materials UJ\. The first one is a 
single barrier (SB) with barrier width 6 = 8 nm and bar- 
rier height V — 0.23 eV. The second system is a double- 
barrier resonant tunneling structure (DB) with barrier 
width b = 5.0 nm, well width w — 5.0 nm and bar- 
rier heights V = 0.23 eV. The third system refers to 
a quadruple-barrier resonant tunneling structure (QB), 
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FIG. 1: (Color on line) Distribution of the first poles of the 
outgoing Green function on the kL plane, SB (full squares), DB 
(full stars) and QB (empty circles), where the parameter L is the 
length corresponding to each system. See text. 



with external barrier widths bi — — 3.0 nm, in- 
ternal barrier widths b2 — b^ — 5.0 nm, well widths 
= i/;2 = W3 = 3.0 nm and barrier heights V = 0.23 eV. 
In all three systems the effective electron mass is taken 
as m = 0.067 me where me is the electron mass. 

For a given potential profile, the parameters of the 
system determine the values of the complex poles 
which are the relevant ingredients to calculate the reso- 
nance states Un{x) and hence the residues r„ appearing 
in both, Eq. ([3]) for the transmission amplitude, and 
Eq. ([50)1 . for the transmitted time-dependent solution. 
Although the procedure to calculate the complex poles 
is known, for completeness, we present in Appendix B, 
a procedure to obtain the necessary number of complex 
poles involving the Newton-Raphson method [33]. The 
set of resonance states {un{x\\ may be obtained using 
the transfer matrix method [l| with the outgoing bound- 
ary conditions given by Eq. (|5)). 

It is of interest to stress that a given potential pro- 
file provides a unique set of resonance poles and 
residues that are calculated only once to evaluate 
Eq. (l30)) . This implies that calculations are much less 
time demanding than calculations involving numerical 
integration of the solution given by Eq. ([T]) where one 
has to perform an integration over k at each instant of 
time, particularly if one is interested, as in the present 
work, to evaluate the above solution at very long times 
and distances. 



A. Complex poles and Transmission coefficient 

Figure [1] exhibits the distribution of the first complex 
poles for the SB (full squares) , DB (full stars) and the QB 
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FIG. 2: (Color on line) Shows the exact transmission (full line) in 
comparison with the analytical formula (Eq. [Sj (dot line) for (a) 
The SB system, using an approximation of 300 poles of the outgoing 
Green function; (b) The DB system, using an approximation of 
1000 poles, and (c) The QB system, using an approximation of 4000 
poles. The inset shows the first isolated resonance and miniband 
corresponding at DB and QB systems respectively. 



(empty circles) systems with parameters as given above. 
In order to facilitate a comparison among the distinct 
distributions, the complex poles of each system are mul- 
tiplied by the corresponding total length L, i.e., respec- 
tively, for the SB, DB and QB systems: L = 8.0 nm, 
L — 15.0 nm and L = 25.0 nm. 

Figures dja), mb), and[2l^c), show respectively, for the 



SB, DB and QB systems considered above, a plot of the 
transmission coefficient T(£') — |t(£')p as function of the 
energy E in units of the barrier height V. Each figure 
presents a comparison between an exact numerical calcu- 
lation using the transfer matrix method (full line), and 
that obtained using the resonance expansion given by Eq. 
^ (dotted line). One observes that the calculations are 
indistinguishable from each other if one considers the ap- 
propriate number of poles, as indicated for each case in 
the caption to Fig. [51 

Let us comment briefly some features of each of the 
above figures. 

Figure[2Ua) exhibits a broad overlapping resonance just 
above the barrier height that is related to the presence 
of the complex energy pole, Ei = £i — iTi/2, with values 
£i = 0.2885 eV and Fi = 0.1045 eV. 

Notice in Figure [IJb), describing a DB system, the 
existence of a sharp isolated resonance in the tunneling 
region with a typical Breit-Wigner or Lorentzian shape 
as exhibited by the inset. The corresponding resonance 
energy parameters are: £i = 0.08 eV and Fi = 1.0278 
meV. At energies above the barrier height the DB sys- 
tem exhibits some transmission resonance structures that 
tend to disappear as the energy increases. 

Figure[5Kc), involving the QB system, exhibits a triplet 
of overlapping resonances along the tunneling region, dis- 
played enlarged in the inset. The corresponding reso- 
nance energies are, respectively, £i = 0.1199 eV, £2 = 
0.1309 eV and £3 = 0.1450 eV and the correspond- 
ing widths, Fi = 4.6270 meV, F2 = 11.9652 meV and 
F3 = 8.4472 meV. Similarly, as in the case of the DB 
system, the QB system exhibits transmission resonances 
above the barrier height as the energy increases. Notice 
that the triplet of overlapping resonances corresponds to 
the first triplet of resonance poles exhibited in Fig. [1] 
(empty circles). This triplet of resonance poles suffices 
to reproduce the transmission coefficient around the cor- 
responding energy range j3Q] . 

B. Time evolution of the transmitted probability 
density 

Let us now investigate the time evolution of the trans- 
mitted probability density \4>{x,t)\'^ using Eq. ([50]) as 
time evolves for different values of a; = Xd. We find 
convenient to plot the dimensionless quantity p{x, t) ~ 
(7\ip{x,t)\'^ in units of t/r, where r stands for the longest 
lifetime of the system, i.e., t = h/Tmin, with Tmin the 
smallest energy width. 

The parameters of the initial cutoff Gaussian wave 
packet, defined by Eq. ((T2|), are 

Xc — —5.0mii, (T = 0.5nm. (38) 

These values give |a;c|/(2cr) ~ 5.0, which implies that the 
condition given by Eq. (fT5)l is satisfied, and hence the 
applicability of Eq. ([50)1 , to calculate the time evolution 
of the transmitted probability density. Notice that a < 
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FIG. 3: (Color on line) Probability density as function of the time, 
in units of Tgg = h/T\, for the SB system (full line). The energy of 
the initial Gaussian state is Eq = V/2. As a comparison, the free 
Gaussian evolution (dotted line) is also plotted. The calculation 
is performed at (a) = 2L short, (b) = 20L medium and 
(c) a;£j = 2 X 10^ L long, distances, where L is the length of the 
SB system. The inset shows the calculation in a semi-ln scale. 
The exact calculation by numerical integration using Eq. ((TJ is 
displayed also in (a) (dashed line). See text. 



L for all the systems considered, i.e., Lsb = 8.0 nm, 
Ldb ~ 15.0 nm and Lqb = 25.0 nm. Also, we choose 
Eq = V/2 for the SB, Eq = £i for the DB system and 
Eq = f 2 for the QB system. The values of the natural 
time scale r are tsb — 5.57 fs for the SB system, tdb — 
0.64 ps for the DB system, and tqb — 0.14 ps for the QB 
system. 

Figure [3] shows a plot of p{xd,t) vs t/r for the SB 
system (full line) for different values of Xd that corre- 
spond to short, medium and long distances from the in- 
teraction region, (a) Xd — 2L, (b) Xd — 2QL, and (c) 
= 2 X lO^L. The results are compared with the cor- 
responding free propagation of the cutoff Gaussian wave 
packet (dotted line). The distinct distances represent dif- 
ferent time scales in the time evolution of the transmitted 
wave packet. Notice that in all the three cases at short 
times, the profiles of the free and transmitted wave pack- 
ets are essentially the same. This follows by noticing that 
the large over-the-barrier energy components of the wave 
packet that impinge on the potential barrier are transmit- 
ted without suffering an appreciable change, as exhibited 
by the behavior of the corresponding transmission coef- 
ficient displayed in Fig. [H (a) {T{E) 1 as E/V > 1), 
and by using a Fourier transform argument that indicates 
that short times correspond to large energies. Figurc[3I^a) 
shows that after reaching its maximum value the trans- 
mitted wave packet (solid line) decays faster than the 
free evolving wave packet (dotted line) . As shown by the 
inset to Fig. ^a), this is so because that time span is 
dominated by the exponential decay of the first top res- 
onance, which in fact after a number of lifetimes suffers 
a transition into a nonexponential behavior as an inverse 
t^'^ power of time as follows from Eq. ([55)) . Figures 
[3jb) and [Sljc) exhibit the time evolution of the proba- 
bility density, respectively, at medium, Xd = 200L, and 
large, Xd = 2 x lO^L, distances and hence medium and 
long times. The corresponding probability density pro- 
files are very similar in both figures. The corresponding 
insets are also similar and show that the time evolution 
goes as the inverse power . At larger distances the 
profile exhibited by Fig. [3] (c) remains unchanged. One 
sees that the profile reflects the energy spectra of the sys- 
tem as shown by a comparison with Fig. [2] (a) for the 
transmission coefficient vs energy. 

Figure m exhibits an analogous calculation of the trans- 
mitted probability density for the DB system (full line) 
and its comparison with the free evolving wave packet 
(dotted line). In this case, at short distances, Xd = 2L, 
the inset shows that the profile of the transmitted wave 
packet is dominated by the transition from the first top 
resonance (See Fig. [UJb)) into the sharp isolated reso- 
nance seated inside the system. At medium distances, 
Xd ~ 200i, one observes a small peak structure around 
t/TDB ~ 8. As the inset displays, it corresponds to 
the exponential decay of the sharp isolated resonance 
situated inside the system. This situation is similar to 
that discussed by Wulf and Skalozub, who considered 
the propagation of a Gaussian pulse near a resonance 
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FIG. 4: (Color on line) Probability density as function of the 
time, in units of tob = ^i/Ti, for the DB system (full line). The 
energy of the initial Gaussian state corresponds to the resonance 
energy, i.e., Eq = £i. As a comparison, the free Gaussian wave 
packet (dotted line) is plotted. The calculation is performed at (a) 
Xd = 2L short, (b) Xd = 200L medium and (c) Xd = 2 X lO^L 
long, distances. The parameter L stands for the length of the DB 
system. The inset shows the same calculation in a semi-ln scale. 
The exact calculation by numerical integration using Eq. ([TJ is 
displayed also in (a) (dashed line). See text. 




FIG. 5: (Color on line) Probability density as function of the time, 
in units of the lifetime Tgg = h/Fi, for the QB system (full line). 
The energy of the initial Gaussian state is equal to the resonance 
energy Eq = £2 ■ As a comparison the free Gaussian wave packet is 
plotted (dotted line). The calculation is performed at (a) Xd = 2L 
short, (b) Xd = 200L medium and (c) ai^j = 2 X 10^ L long distances, 
where L is the length of the QB system. The inset shows similar 
calculations in semi-ln scale. The exact calculation by numerical 
integration using Eq. JTJ is displayed also in (a) (dashed line). 
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level . The inset shows that eventually at longer times 
there is a transition to nonexponential decay as an in- 
verse power of time. Finally at very large distances, 
Xd = 2x 10^ L and very long times, of the order of IO'^td b , 
in a similar fashion as in the previous system, the pro- 
file of the transmitted wave packet reflects already the 
structure of the energy spectra of the DB system. The 
corresponding inset to Fig. UJ^c) shows that the sharp 
structure around I/tdb ~ 7.2 x 10^ evolves at long times 
in a nonexponential fashion. 

Figure [5] exhibits also a situation similar to the exam- 
ples discussed above, for the time evolution of the trans- 
mitted probability density of the QB system. Again Figs. 
[5] (a), (b) and (c) refer, respectively, to short, Xd — ^L, 
medium Xd — 200L, and large, Xd — 2 x lO^L, distances. 
At short distances, Xd = 2L, it is worthwhile to notice 
the presence of Rabi oscillations in a similar fashion as 
occur in the decay of multibarrier systems [ssj. These 
oscillations represent transitions among the closely lying 
resonance levels of the QB system. Again as the distance 
and the time increase, the resonance levels decay, first 
exponentially and then nonexponentially, as depicted in 
the inset to Fig. Efb). At still larger distances the decay 
is purely nonexponential, as an inverse cubic power of 
time, as shown by the inset to Fig. \^c). Notice that in 
Fig. [njc) the profile of the transmitted wave packet re- 
sembles already the energy structure of the transmission 
coefficient. 

It is of interest to compare our results with the case 
of a cutoff incident plane wave impinging on a multi- 
barrier system. This case, corresponding to the limit 
of an infinitely broad Gaussian wave packet, has been 
considered recently by Villavicencio and Romo (36j . us- 
ing the formalism developed in Ref. 3, to investigate 
the propagation of transmitted quantum waves in these 
systems. There, for incidence energies Eq below the low- 
est resonance energy of the multibarrier system, a se- 
ries of propagating pulses (forerunners) are observed in 
the transmitted solution traveling faster than the main 
wavefront. It is shown that each forerunner propagates 
with speed v{£n) =[2m£„/m]^/^ associated with the nth 
resonance of the system, thus establishing a relation- 
ship between the sequence of forerunners and the reso- 
nance spectrum of the system. However at asymptot- 
ically long times the forerunners fade away, since the 
solution ■tp{x,t) ~ t(fco) exp(ifcoa;) exp(— ii5oi)/^)j with 
fco — [2mEo]^^'^ /h. This yields for the transmitted prob- 
ability density i)P = K(fco)P, a result very different 
from the case of Gaussian wave packets of finite width 
considered here. 
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C. Reconstruction of the energy spectra 

In order to exhibit more clearly the relationship be- 
tween the time evolution of the transmitted Gaussian 
wave packet and the energy spectra of the system, 
pointed out in the previous subsection, one may proceed 



FIG. 6: (Color on line) Figures (a), (b), and (c) refer, respectively, 
to the single barrier (SB), double-barrier (DB) and quadruple- 
barrier (QB) systems. In each figure: T{E) stands for the transmis- 
sion coefficient in units of Eq , the energy of the incident wavepaket 
(solid line); (^(x, to) defined by Eq. 1391 1 in units of the parameter r], 
for distinct values of to , as specified in the insets to each figure. As 
time increases the transients end by reproducing the energy spectra 
of the corresponding systems. See text. 
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as follows. First, instead of evaluating the transmitted 
probability density by fixing x = Xd and varying the time 
t, i.e., \ij{xd, as discussed in the previous subsection, 
we consider instead a fixed value of the time t ~ to and 
vary x, i.e., \ip{x,to)\'^ . It is not difficult to see that 
the plot of \ip{x,to)\'^ vs X looks identical to the specu- 
lar image, with respect to the vertical axis at the origin, 
of \'4'{xd,t)\'^ vs t. Second, in analogy with the calcula- 
tion of the transmission coefficient in the energy domain, 
we divide \'ip{x,tQ)\'^ by the free evolving Gaussian wave 
packet \'4'[ix,to)\'^, given by Eq. (^0]) . We define the 
quantity (^{x, io) as the ratio of these quantities, namely. 



(39) 



Third, it is convenient to plot the transmission coefficient 
T{E) in units of E/Eq, with Eq, the incident energy of 
the corresponding Gaussian wave packet. This allows to 
relate the values of E/Eo with the values of a parameter 
77 defined as 



1 2 



xq- L 



(40) 



as follows. The above expression for rj is based on the 
argument that for x ^ L, E = U^k'^ /2m with hk/m = 
{x—L)/to, and i?o = ?i^fcQ/2TO with hko/m — {xQ — L)/tQ, 
where xq — L in the distance that a free particle travels 
in time to. Hence, E = riE^. Figures [6] (a), (b) and 
(c) exhibit respectively, for the SB, DB, and QB sys- 
tems, considered in the previous subsection, the plots of 
C(a;,to) vs rj. For each of the above figures, three graphs 
are plotted. As indicated in the inset to each figure, each 
graph of C.{x, to) corresponds to a distinct value of to and 
hence oi xo — L. The above figures also exhibit a plot of 
T{E) in units of E/Eo (solid line). Notice, as pointed out 
previously, that the value of Eo differs for each system. 
One may appreciate, in each figure, the transient behav- 
ior of the transmitted Gaussian wave packet. For small 
values of to, (,{^^^0) just reproduces the fastest compo- 
nents of the energy spectra of the corresponding system 
and as to increases it goes into a transient behavior that 
ends when the transmission energy spectra of the system 
is reconstructed, as shown analytically by Eq. p7|) . 



D. Remark on the tunneling time problem 

Our results are of relevance for the tunneling time 
problem [ll[. Here, the question posed is: How long 
it takes to a particle to traverse a classical forbidden 
region? One of the approaches considered involves the 
tunneling of wave packets. Here one usually compares 
some feature of the incident free evolving wave packet 
(usually a Gaussian wave packet) and a comparable fea- 
ture of the transmitted wave packet, commonly the peak 
or the centroid, and a delay is calculated. Many years 



ago, Biittiker and Landauer [37| argued that such a pro- 
cedure seems to have little physical justification because 
an incoming peak or centroid does not, in any obvious 
causative sense, turn into an outgoing peak or centroid, 
particulary in the case of strong deformation of the trans- 
mitted wave packet. Our results for the transient behav- 
ior of the transmitted wave packet supports that view in- 
dependently of whether or not there is initially a strong 
deformation of the transmitted wave packet. Even if the 
transmitted wave packet is initially no deformed, as in 
Fig. El (a) for a single barrier system, as time evolves the 
profile of the transmitted Gaussian wave packet varies to 
finally reproduce the energy spectra of the system and 
hence there is no a unique way to answer the question 
of how long it took to the initial packet to traverse the 
system. 



V. CONCLUDING REMARKS 

The main result of this work is given by Eq. ([HU)) . 
which provides an analytical solution to the time evolu- 
tion of a Gaussian wave packet along the transmission 
region for scattering by a finite range potential in one 
dimension. We have focused our investigation to cases 
where the Gaussian wave packet is initially far from the 
interaction region, i.e. fulfills Eq. (US]), and is suffi- 
ciently broad in momentum space so that all sharp and 
broad resonances of the system are included in the dy- 
namical description. We have obtained analytically and 
exemplified numerically for single and multibarrier quan- 
tum systems, that the profile of the transmitted Gaus- 
sian wave packet, exhibits a transient behavior that at 
large distances and long times becomes proportional to 
the transmission amplitude of the system, i.e., Eq. p7p . 
This predicts the final destiny of the transmitted wave 
packet in a coherent process. It is also worth to empha- 
size that the analytical expression for the transmitted 
wave packet yields, at a fixed distance and asymptoti- 
cally long times, a t~^/^ behavior with time, i.e., Eq. 
psp . This result corroborates numerical calculations for 
Gaussian wave packets colliding with square barriers and 
extends previous analysis to arbitrary potentials of finite 
range One should recall that the set of poles {k„} 
and residues {r„}, which is unique for a given poten- 
tial profile, is evaluated just once to calculate the time- 
dependent solution given by Eq. ((30l) . The number of 
poles required for a dynamical calculation corresponds 
to the number of poles necessary to reproduce the exact 
transmission amplitude using Eq. This is in con- 

trast with calculations involving numerical integration of 
the solution, using Eq. ([T]), where one has to perform 
an integration over k at each instant of time and hence 
the calculation is much more demanding computationally 
particularly at large distances and long times. Further 
work is required to extend the results of the present in- 
vestigation to wave packet dynamics in multidimensional 
tunneling [39j . Our analytical solution for the transmit- 
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ted wave packet might be of interest in connection with 
the long debated tunnehng time problem. 



Then, the nonexponential contribution of each pole k„ 
in Eq. ([XT]) reads 
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APPENDIX A: ANALYSIS OF Lu{-iz). 

Here we show that the contribution of the term a;(— iz), 
appearing on the right-hand side of Eq. p7|) , to the time 
evolution of the transmitted Gaussian wave packet may 
be neglected provided the condition given by Eq. (fT8|) 
is fulfilled. The contribution corresponding to Lu{~iz) 
reads, 



(27r)V4^2. 



2tt 



(Al) 

where z, defined by Eq. (jlSp . is written as z = i(k — 
kQ)a with kg = ko — ixc/2a'^. In general, it is necessary 
to calculate numerically the integral term given by Eq. 
(lAip . However, for the particular case specified by Eq. 
(fT5)l . i.e., \xc/2a\ > 1, that implies that \z\ > 1 for aU 
values of k, one may use the asymptotic expansion of the 
Faddeyeva function u!{—iz) [3ll |33| 



uj{—iz) 



I ^ r(j + 1/2) 

^jr:,[{k-k-)aY=+^ 



N 



TT ^ 

l2j 



(2j)!cr2i+i Kk- 



k' ' 



(A2) 



where the quantities D,, denote a 27-th derivative oper- 
ator. 

Substitution of Eq. (IA2p into Eq. (jAip allows to ex- 
press each integral term in the sum as 



2^ 



dfc- 



(k k^^^k Kji) 
where we have used the identity 
1 1 



k' 

ftp 



(fc kQ)[k Kn) 



k kr) k 



(A3) 



(A4) 



and the arguments of the Moshinsky functions M{y'o) and 
M{yn) are given respectively by. 



-jTr/4 



and 



fik'^ 



(A5) 



(A6) 



{2^Y'^^ 1 A r(j + 1/2) 



2 3 + 1 



Af(y^) - M(y„) 



(A7) 



The dominant term in powers of a in Eq. (jA7p occurs 
for J = and hence the nonexponential contribution of 
each pole is given by 



2V'^ 1 Af (y(,) - Af(2/„) 



(A8) 



Recalling that the factor ^JLo{^ZQ) = exp(zQ /2)^J erfc(zo) 
and that zq ^ — 1 (Sll [s^ one obtains. 



(A9) 



It follows then, by substitution of (jA9p into (jASp and 
comparing the resulting expression with Eq. (I23|) , taking 
into account that the corresponding Moshinsky functions 
yield contributions of the same order of magnitude, that 



\rn{x,t)\ 



\rn{x,t)\. 



(AlO) 



The above expression demonstrates that provided Eq. 
(flS)) is satisfied, the nonexponential contribution i/jJJ'^ 
may be neglected. 



APPENDIX B: CALCULATION OF COMPLEX 
POLES OF THE TRANSMISSION AMPLITUDE. 

It is well known that the transmission amplitude t(/c) 
for a potential V(x) of finite range, i.e., extending from 
X — Q to X — L, possesses an infinite number of complex 
poles K„ that in general are simple [28l |. These complex 
poles correspond to the zeros of the element t22{k) of the 
corresponding transfer matrix 



t(fc) 



i22(fc) 



(Bl) 



The set of complex poles of t(fc) may be calculated 
using the Newton-Raphson method [33j. This method 
approximates a complex pole /t„ by using the iterative 
formula 



^22«) 
^22 



(B2) 



where i22(^) — dt22(k)/dk. The approximate pole kJ^"*"^ 
goes into the exact pole, at a given degree of accuracy, as 
the number of iterations increases. In order to apply this 
method, it is necessary to provide an appropriate initial 
value for the approximate pole k^. 
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FIG. 7: (color on line) (a) Comparison of the transmission coeffi- 
cient T{E) as a function of energy in units of the potential height 
Vb for the QB system of the exact calculation using the transfer 
matrix method (solid line) with resonance calculations using Eq. 
S3l for a number of poles: 10 (dashed line), 100 (dash-dot line) and 
1000 (dotted line), (b) A similar calculation for the transmission 
amplitude Ret(fc)t)sImt(A;). 



In general for systems formed by a few alternating bar- 
riers and wells, as exemplified by Fig. [21 the transmis- 
sion coefficient vs energy may be roughly characterized 
by three regimes: Regime I, characterized by sharp iso- 
lated resonances (as in Fig. [U^b)) or groups of well de- 
fined overlapping resonances (as the resonance triplet in 
Fig. (He)). This regime occurs usually for energies below 
the potential barrier height and refers to complex poles 
that are seated close to the real fc-axis; Regime II, char- 
acterized by broad overlapping resonances. This regime 
is commonly found close to the potential barrier height 
and may extend up to energies 3 or 4 times the potential 
barrier height, as exemplified in all Figs. [2 and Regime 
III, involving much higher energies, well above the bar- 
rier height. There the transmission coefficient does not 
exhibit any appreciable resonance structure and just fluc- 
tuates very closely around unity. 

There is in general no analytical expression for any 
initial approximate pole k^. An exception occurs along 
the regime III, where there exists an asymptotic formula 
for the location of complex poles which is valid for very 



mr .2 
--z-ln(n) 



0(1); n>l. 



(B3) 



One may substitute Eq. (|B3p into Eq. (|B2p to obtain 
the pole Kn = an — i/3n for that very large value of n, 
say for example, n = 4000. Equation (|B3|) provides a 
relationship between the real parts of the {n + l)th and 
(n — l)th poles with the nth pole 



Onil « Q!n ± — = a„±i, 

and for the corresponding imaginary parts. 

Hence one may write 

where the step is given by 



(B4) 



(B5) 



(B6) 



(B7) 



Then, one may calculate the (n — l)th pole by substitut- 
ing Eq. (jB6p into the iterative Newton-Raphson formula 
to evaluate the pole k„_i. Repeating this procedure suc- 
cessively allows to generate the poles for smaller values 
of n. Clearly this procedure permits also to obtain the 
poles for larger values of n. As the value of n dimin- 
ishes, however, on may reach a situation where, even if n 
is still large, the iterative Newton-Raphson formula may 
fail. We have found that in this circumstance Eq. (|B7p 
still holds but Eq. (jBSP becomes inaccurate. In order 
to circumvent this situation one may proceed as follows. 
Once, as indicated above, that it is determined that the 
pole K„ is asymptotic and has been calculated, one de- 
fines a rectangular region on the complex k plane 
whose center contains the pole This region is char- 

acterized by 



/„_i = [a„_i - A,./2, a„_i + Ar/2] x 
[-/?„- A,/2,-/3„ + A,/2], 



(B8) 



where Aj is a controllable parameter. Since the imagi- 
nary values of neighboring poles do not differ substan- 
tially, it is sufficient to choose 



(B9) 



If, as indicated above, the iterative formula given by Eq. 
(jB2p fails for a given initial value k'^-i, then a new ini- 
tial value K^-i generated randomly according to the 
expression 



(BIO) 



where, the parameters "/r and 7^ are random numbers 
that vary, respectively, along the intervals —0.5 < jr < 
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0.5 and —0.5 < 7i < 0.5 to guarantee that the gener- 
ated pole hes within the region If the condition 
|^22(k°_i)| < 1 is fulfilled, then the iterative formula (|B2p 
is applied. Otherwise or if the calculated pole lies out- 
side In-i, that pole is disregarded and a new initial pole 
is generated according to the above procedure. Usually, 
after a few random attempts convergence to a new pole 
is obtained. If after many random attempts (M=1000 
for the examples considered in this work) no convergence 
is achieved, that may suggest that Regime II has been 
reached. This means that Eq. (|B7[) does not hold any- 
more. Then, it is convenient to define from that pole 
inwards thinner rectangular regions /„_i. For the ex- 
amples considered in this work, we choose Ar = ■n/2QL 
and for Ai — 2/3„. Clearly, in this case some rectan- 
gular regions do not possess any poles. This procedure 
is capable to generate also the poles in Regime I. Al- 
though in Regimes I and // the above procedure may 
generate repeated poles, a consequence that Eq. (|B4p 
does not hold, these poles may be easily identified and 
disregarded. For Regime I there is the alternative simple 
procedure to generate the initial values by the rule 
of the half-width at half-maximum of the Breit-Wigner 
formula for the transmission coefficient. 

Once a set of TV complex poles {k„} has been obtained, 
one may evaluate the transmission amplitude given by 
Eq. ([2]), by running it from —N to N. One might 
then make a comparison of the resonance expansion, for 
different values of the number of poles, with the exact 
numerical calculation using the transfer matrix method 



m to establish the appropriate number of poles for a 
given energy interval. Figure Elja) provides a plot of the 
transmission coefficient vs energy for the QB system dis- 
cussed in the text for the exact numerical calculation 
using the transfer matrix method (solid line) and reso- 
nance expansions of t(fc) for distinct number of poles: 
iV 10 (dashed line), N = 100 (dash-dot fine) and 
N — 1000 (dotted line). The energy interval extends up 
to 5 times above the barrier height and one sees that as 
the number of poles increases the agreement with the ex- 
act calculation becomes better. Notice that already with 
N ~ 100 poles, the transmission coefficient is well re- 
produced for energies below the potential barrier height. 
Notice also that the calculation involving 1000 poles is 
still slightly different from the exact calculation in the 
interval 4.0 < E/Vq < 5.0. The calculation for the same 
system presented in Fig. [2tc), that involves 4000 poles, 
is indistinguishable from the exact calculation. One sees 
that away from sharp resonances, more resonance terms 
are required to reproduce the exact calculation. This is 
particularly striking in energy intervals where T{E) fluc- 
tuates very close to unity where a very large number of 
resonance terms is necessary to reproduce the exact cal- 
culation. Fortunately, very distant resonance poles are 
not difficult to calculate. Figure [Tl^b) exhibits similar 
calculations for the transmission amplitude. Here it is 
plotted Ret(fc)wslmt(fc), to show that the resonance ex- 
pansions of the transmission amplitude become closer to 
the exact calculation as number of poles in the calcula- 
tion increases. 
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